Coexistence of the superconducting energy gap and pseudogap above and below the 
transition temperature of superconducting cuprates 
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We express the superconducting gap, A(T), in terms of thermodynamic functions in both s- and 
d-wave symmetries. Applying to Bi2Sr2CaCu2 08+5 and Yo.8Cao.2Ba2Cu307_i we find that for all 
dopings A(T) persists, as a partial gap, high above Tc due to strong superconducting fluctuations. 
Therefore in general two gaps are present above Tc, the superconducting gap and the pseudogap, 
effectively reconciling two highly polarized views concerning pseudogap physics. 



On cooling a superconductor (SC) below Tc coherent 
pairing of electrons opens a gap. A, centered at the Fermi 
level. In a conventional SC A(T) closes at Tc but for un- 
derdoped cuprates a partial gap is found to persist above 
Tc and this is widely attributed to the so-called pseudo- 
gap JL]. The field is sharply divided as to the origins of 
the pseudogap. One view is that it is some form of pre- 
cursor SC state while another is that it arises from some 
correlation that competes with the SC state [if, so that 
the two gaps coexist below Tc. The inherent physics for 
each scenario is fundamentally different. In the former 
case a phase-incoherent SC state [i| emerging from RVB 
physics high above Tc Q is often invoked, implying a very 
large SC energy gap which falls rapidly with increasing 
doping. In the latter case, it is the pseudogap, arising 
from some independent competing correlation, that has 
the large energy scale and the pseudogap closes abruptly 
at a putative ground-state quantum critical point lying 
within the SC dome at pcrit^O. 19 holes/Cu [j]. 

Because these two scenarios differ so radically it re- 
mains a central challenge to identify the nature of these 
energy gaps. Is there, indeed, one or two distinct gaps? 
Here we present a new method to calculate A(T) from the 
electronic specific heat. We show that, for the cuprates 
at any doping, A appears to remain finite above Tc re- 
flecting a partial gap arising from strong SC fluctuations 
and, in the underdoped region, coexists there with the 
pseudogap. Thus, in a sense, both scenarios prove to he 
correct. There are two gaps above Tc just as there are two 
gaps below Tc so that both fluctuations and competing 
pseudogap correlations play a key role in HTS physics. 

Using a high-resolution differential technique Loram 
et al. p] have been able to isolate the electronic spe- 
ciflc heat from the much larger phonon term in a number 
of high- Tc cuprates. This has allowed many important 
conclusions to be drawn including the fact that, due 
to strong SC fluctuations, the mean-fleld (MF) transi- 
tion temperature, T™^, determined from entropy conser- 
vation, lies well above the observed Tc value (by up to 
50K) 0,13. Hereafter, we drop the descriptor 'electronic' 
and by the terms specific heat, Cp, specific heat coeffi- 
cient, 7 = Cp/T, entropy, S, internal energy, U , and free 



energy, F, we mean the electronic components of these. 

We draw largely on Ferrell [3] and extend to include 
d-wave SC. Starting from the BCS Hamiltonian he shows: 



dTcJj, 



-Co" N{0)TcQ{t) 



(1) 



where 7V(0) is the DOS at the Fermi level, t = T/Tc, 
Q{t) = (A(T)/Ao)^ and we include the additional factor 
C = 1. For an anisotropic gap we take A to be the 
amplitude of the k-dependent gap. In this case Ferrell's 
A{T)^ should be replaced by a Fermi surface average 
(Ak(T)2) = CA(T)2 where C = 1 for s-wave and 1/2 for 
d-wave. The BCS gap ratio a = Aq/ kpTc = (7r/7£;)e'-^ 
where je — 1.781... is Euler's constant and C = for 
s-wave while C — In 2 — 1/2 for d-wave Q . Ferrell then 
integrates Eq. [1] over all Tc > T to effectively obtain: 

AF(T) = T„(T)-T,(T,Tc) 

= CN{0)Al t'-^ Q{t') At' (2) 

Ferrell's intention was to adopt a model T dependence 
of Q{t) from which to calculate AF{T). Our task is the 
opposite, to calculate A(T) from AF{T) derived from 
specific-heat data. By differentiating each side of Eq. [2] 
with respect to T and rearranging we obtain 

CiV(0)A(T)2 = 2AT(T) + TA5(T) 

= 2AC/(T) - TA5(T) (3) 

which expresses A(T) directly in terms of thermody- 
namic functions AT, AS = Sn — Sg and AU = Un — Ug. 

Quite generally, for a second-order MF phase transition 
near Tc, AT(T) = -\A-fc{T - Tcf , so 



A{Tf 



'2T^Ajc 
CN{0) 



(1 - T/Tc) 



(4) 



where A7C is the jump in 7 at Tc. This means that 
the coherence length ^(T) = hVp /TrA{T) has the correct 
(1 — dependence near Tc. 

We have computed (2AC/(T) - TAS'(T))/2Af/(0) for 
both s- and d-wave weak-coupling BCS and in Fig. [I{a) 
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FIG. 1: (Color online) (a) The temperature dependence of the 
square of the normalized SC gap function (A(r)/Ao)'^ for s- 
wave and d-wave weak-coupling BCS (dashed curves). These 
are compared with the values of this parameter calculated 
from Eq,|3] using (2A[/(T) - TA5'(r))/2AI7(0). The inset 
shows the T-dependence of each contribution 2AU and TAS. 
(b) The T-dependence of the SC gap function A(r) calculated 
using Eq. [3]for d-wave strong-coupling BCS with 2a = 4.28 
(weak coupling), 5, 6 and 7. 



we compare these (solid curves) vifith the theoretical T- 
dependence of (A(T)/Ao)^ (dashed curves). For both 
symmetries there is excellent agreement across the entire 
T-rangc and the gap amplitude satisfies: 



A^^(O) = A;7(0) = C/o - -C^(0)A2. 



(5) 



This is just the ground-state condensation energy. The 
inset shows the individual contributions 2AU and TAS 
to A(T). AU{T) passes through a maximum while sub- 
traction of the entropy term recovers the canonical mono- 
tonic T-dependence of the s- or d-wave gap. 

Ferrell's theory is strictly for weak-coupling BCS, 
based on the logarithmic relation between the pairing in- 
teraction and Tc. Extending to strong-coupling we may 
employ the Padamsee a-model approximation [10] where 
the ratio a = Ao/fcsTc is the only adjustable parame- 
ter and A(T)/Ao is assumed to follow the weak-coupling 
BCS form for all a. As a cancels in Eq. [2] we might 
still consider using Eq. [3]to calculate A(T). In the case 
of Pb, a strong-coupling superconductor, we have calcu- 
lated AT" and AS" from critical-field measurements, and 
thence A(T) using Eq. [H We find excellent agreement 
with measurements of the gap from tunneling including 
the flattening of A(T) relative to the BCS T-dependence. 
This gives us confidence to extend beyond weak-coupling. 



as may be necessary for the cuprates. 

Accordingly, we used the a-model to calculate AF and 
AS* for 2a = 4.28 (weak coupling), 5, 6 and 7 in a d-wave 
scenario, employing the same method as Padamsee et al. 

Fig.IUb) shows A(T) calculated for each case usmg 
Eq. [3l The fine black curve under the blue dashed curve 
for 2a = 4.28 is the theoretical weak-coupling BCS gap, 
for which the match is exact. In the figure 7„ is the 
normal-state (NS) value of 7, which is assumed to be 
T- independent. In strong coupling, 7^ is enhanced by a 
factor (1 -t- A) above its Sommerfeld value, viz. 



In 



^2fc|iVo(l-HA), 



(6) 



where A is the usual electron-boson coupling parame- 
ter in Eliashberg theory ll| and iVo is the bare band 



DOS, un-renormalized by electron-boson or Coulomb ef- 
fects. Thus A(T) in Fig. [Ijb) is expressed in units of 

Y^^TrfesTc-v/l + A. Leaving aside the absolute magnitude 
of A, the T-dependence of A evidently fiattens with in- 
creasing coupling. Though this violates the main premise 
of the a-model, the a-model could potentially be refined 
by calculating A(T) iteratively. For the time being, Eq.[3] 
seems to be a satisfactory approximation for strong cou- 
pling if 2a<5, as we find for the cuprates [l^j . 

We now apply this analysis to the electronic specific 
heat of the HTS cuprates reported by Loram et al. [Hj. 
Fig.EJa) shows the previously reported analysis Q of the 
specific heat coefficient, 7(T), for Y0.8Ca0.2Ba2Cu3O6.75 
used to determine the mean- field Tc value, T^^ . At 
this doping {p — 0.185) the pseudogap is absent and 
the NS coefficient jn{T) is essentially constant (dashed 
line), the MF 7 in the SC state deduced by 

entropy balance, namely the area abc equals the area 
cde. Also by entropy balance the grey shaded area un- 
der the fiuctuation contribution equals the hatched area 
which therefore defines T"^^ . By integrating 7(T)— 7„(T) 
in Fig. ma) we obtain AS" = Ss — Sn and similarly 
A5""/ = /o^(7™/ -7„)dT. These are plotted in Fig.lHc) 
by the solid and dashed curves, respectively, where only 
every 4''' data point is shown. These may in turn be inte- 
grated to generate AF{T) and AF™^ (T) and these com- 
bined with TAS{T) and TAS""/(T) to generate A(T) 
and A'"^f[T) using Eq. El These are plotted in Fig. ^h) 
where -/V(O) is obtained using Eq. [HI with A = 0. The 
actual gap should be larger by the factor \/l + A. 

The first point to note is that A"^f (T) is found to fol- 
low almost precisely the BCS temperature dependence. 
This means that the HTS systems are close to weak- 
coupling behavior as we have previously deduced Q thus 
justifying the basic assumptions of our analysis. Even if 
A is appreciable one could invoke the Padamsee approach 
to renormalize the magnitude of A provided 2A/kBTJ!^^ 
does not greatly exceed the BCS value of 4.3. Secondly, 
with increasing temperature A(T) starts to fall below 
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FIG. 2: (Color online) (a) Reproduced from analysis of 
the specific heat coefficient, 7(r), for Y0.8Ca0.2Ba2Cu3O6.75 
to determine the mean-field Tc value, T™-^, showing the de- 
duced MF coefficient 7™-'^ and the symmetric fluctuation con- 
tribution (grey shading). By entropy balance the hatched 
area equals the shaded area under the fluctuation term, (b) 
solid curve: the SC energy gap Ao calculated using Eq. [3]and 
dashed curve: its MF value, A™-'^ calculated from 7"'-^ in (a), 
(c) solid curve: the entropy difference AS* = Ss — Sn calcu- 
lated by integrating (7 — 7n) from (a); dashed curve AS""-'^ 
calculated by integrating (7'"-'^ — 7,1) from (a). 



A'"/(T) at the onset of SC fluctuations below T^. At 
there is an inflexion in A(r) which then remains finite 
and falls only slowly to zero above T^. As it does so it 
becomes less well defined due to the square root in Eq. |31 
At Tc the coherent SC state vanishes and this finite resid- 
ual "gap" reflects a fluctuation-induced loss, above Tc, of 
spectral weight in the DOS at Ep - just as described by 
Fig. 10.2 in Larkin and Varlamov [l3|. 

A similar analysis was carried out for many other 
doping levels for both Yo.8Cao.2Ba2Cu307_5 and 
Bi2Sr2CaCu208-i-5. Now we must take into account the 
effects of the proximate van Hove singularity (vHs) on 
the overdoped side (7n(T) rises with decreasing T) and 
the pseudogap on the underdoped side (7n(T) falls with 
decreasing T). To do this we still use entropy balance 
but employ a rigid ARPES-derived dispersion, which im- 
plicitly contains the vHs, to determine the doping evo- 
lution of the background 7n(T) and Sn{T). We use the 
model of Storey et al which includes a non-nodal 
NS pseudogap at (7r,0) reflecting the formation of hole 
pockets as described, for example, by the Fermi-surface- 
reconstruction model of Yang, Rice and Zhang 15| . The 
pseudogap closes abruptly at p « 0.19. For more details 
see SM For Bi2Sr2CaCu208+5 the data for S{T) 
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FIG. 3: (Color online) The T-dependence of (a) the conden- 
sation free energy AF{T) and (b) A obtained using Eq. [S] 
for Bi2Sr2CaCu2 08+5 . Curved arrows show increasing dop- 
ing from p = 0.12 to p = 0.22. Note that the absolute value of 
A(r) is larger than that shown by a factor of vT+A. Inset; 
the doping dependence of the BCS ratio AF{0)/jn{Trff de- 
termined for Yo.8Cao.2Ba2Cu307_5. The BCS value of 0.17 
is preserved for p > 0.19 but falls rapidly with the opening of 
the pseudogap. 



and the dispersion-derived Sn (T) are already reported by 
Storey et al. Q. By integrating AS'(T) = S'(T)-5„(T) 
we obtain AT(T) which is shown in Fig.|3Ua) for 11 dop- 
ings spanning under- and over-doped regions. 

From AF{T,p) and TAS{T,p) we calculate A(T,p) 
using Eq. [3l This is plotted in Fig. ^h) and as before 
the deduced gap does not vanish at Tc. Rather, it inflects 
there and then persists some 20 K above Tc in overdoped 
samples and up to 70 K above Tc for underdoped samples. 
Residual gaps above Tc are perhaps not new however 
they are usually confused with the pseudogap [16]. We 
distinguish between the two gaps, as follows. 

The residual gap that we observe above Tc arises from 
SC fluctuations near Tc which are distinguished by a fluc- 
tuation term in 7(T) which is symmetric over a narrow 
range about Tc (see grey shaded areas in Fig.[2ja)). The 
pseudogap is altogether different. Its effects are not cen- 
tered on Tc but extend over a broad temperature range 
up to 300 K or more, and is distinguished by: 

(i) a broad suppression of S{T)/T as T is reduced 1, 
[l3 | , corresponding precisely to the suppression of the spin 
susceptibility, Xs(T), long observed in NMR [13] ■ 

(ii) the abrupt reduction in the jump A7C at Tc with 
the opening of the pseudogap at p = 0.19 holes/Cu. As p 
is reduced below 0.19 A7C is rapidly diminished, reflect- 
ing a crossover from strong to weak superconductivity. 
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FIG. 4: (Color online) False color plot of (a) the SC gap, A{T,p) and (b) the BCS-normalised condensation free energy 
AF (T) /■y„{T^-^ )^ across the {p — T) phase diagram for Yo.8Cao.2Ba2Cu307_a. In (a) the color scale runs from 20 meV (red) 
to meV (blue) and in (b) from 0.18 (red) to (blue). Also shown is the observed Tc and determined previously [^. 

The SC gap extends well above Tc while AF{T) is cut off at Tc- Also shown is T*(p) (white line) previously reported for 
Yi_a;Caa;Ba2Cu307_5 (epitaxial thin films: down-triangles; poly crystalline: up-triangles) and also as reported by Daou et al. 
[2^ (magenta squares) for YBa2Cu307-5 single crystals. In (a) the dashed white curve is the envelope of the pseudogap in 
the underdoped region and the residual SC gap in the overdoped region. In (b) Af (0)/7„(T<r^)^ adopts the BCS value 0.17 
across the overdoped region but collapses rapidly at T* confirming that the T* line does indeed terminate at p = 0.19. 



(ill) a relative insensitivity to the effect of a magnetic 
field or impurities [3, [l^ in distinct contrast to the pair- 
ing gap arising from SC fluctuations. 

As A(T) persists above Tj,, even in overdoped samples 
where the pseudogap is absent, it must therefore arise 
from SC fluctuations above Tc- On theoretical 13[ and 
experimental [lO] grounds, this will cause a gap-like loss 
of spectral weight and an associated entropy loss which 
underlies the residual A(r). In further support, Gomes 



et al. [2l| observe a spatially inhomogeneous partial gap 



above Tc in tunneling spectroscopy up to a temperature 
Tp.max which closely matches our T^-^ . Such a gap is 
also seen in ARPES [22|. This partial gap also probably 



underlies the anomalous Nernst effect observed in under- 
and over-doped samples between Tc and T^^ [23| . 

Returning to Yo.8Cao.2Ba2Cu307_5, similar results 
are found. Fig. IDJa) shows a false-color plot of the mag- 
nitude of A(r, p) across the p — T phase diagram, along 
with Tc{p) and the previously determined T™f{p) 
A finite gap extends above Tc and indeed above T™^, 
reflecting again the presence of strong SC fluctuations, 
though neither the gap nor T^f extend as high as in the 
case of Bi2Sr2CaCu208+5. To emphasize the crucially 
important distinction between pseudogap and SC corre- 
lations we also plot in Fig. |4ja) previously determined T* 
values for Yi_a;Caa;Ba2Cu307_5 (epitaxial films: down- 
triangles, polycrystalline: up-triangles). T* was deter- 
mined in the usual way by the downturn from linear re- 



sistivity, with the added precaution of using a magnetic 
field to distinguish between SC fluctuations and the pseu- 
dogap near Tc [H, [3l ■ T* cuts through the crescent of 
finite A(r, p) above Tc, falling to zero at p « 0.19. 

This issue continues to be debated. Many groups es- 
pouse a T* line similar to the white dashed curve in 
Fig. Ufa) that extends above the SC dome, across the 
overdoped region. Daou et al. [23 | are a recent exam- 
ple. We therefore also plot Daou's T* data points (ma- 
genta squares) for YBa2Cu307_5 and they are in excel- 
lent agreement with our data shown by the white trian- 
gles and solid white curve. The white dashed curve is the 
envelope above Tc of a finite gap-like feature, whether the 
SC gap or the pseudogap. Unless these gaps are distin- 
guished it is not surprising that many groups have failed 
to see that T* terminates abruptly aX p k, 0.19. 

Definitive evidence for the termination of T*(jj) at 
p = 0.19 is shown in Fig. [H^b). Here we plot a false 
color plot of the ratio AF{T)/^n{T^^Y across the phase 
diagram. Note that we have used T™^ and not Tc as 
the normalising energy scale. The value of this ratio at 
T = is also plotted in the inset to Fig. [SJa). The uni- 
versal BCS c?-wave value for this ratio is 0.17 and, indeed, 
this value is obtained across the entire overdoped region 
for p > 0.19. But with the opening of the pseudogap for 
p < 0.19 the ratio is seen in Fig.lDJb) to collapse abruptly, 
clearly delineating the termination of the pseudogap T* 
line. Moreover, this shows that the pseudogap and su- 
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perconducting gap coexist below Tc as they do above Tc- 
We are thus obhged to conchide that T*{p) cuts the SC 
dome and terminates at p w 0.19, contrary to the infer- 
ence of Daou et al. though in fact their data is fully 
consistent with our scenario. 

Finally, the magnitude of Aq obtained here is a lit- 
tle lower than that observed previously by e.g. infrared 
measurements [m where the amplitude is about 25 meV. 
But recah that A(r) in Figs. [3] and H is yet to be 
enhanced by the factor vT+A. 

In summary, we have shown that the SC gap, A(T), 
may be calculated from the electronic specific heat and 
we apply to the cuprates. For all dopings a residual finite 
A(r) extends up to 70K above T^, reflecting a fluctuation- 
induced loss of spectral weight at Ep. This crescent of 
residual SC gap above Tc is cut by the T* line showing 
that two gap- like features are present above Tc, one ex- 
tending across the entire SC phase diagram due to strong 
SC fluctuations and the other present only in the optimal 
and underdoped region due to the pseudogap. The ratio 
C/o/7«(TJ"/)2 adopts the BCS weak couphng value (0.17) 
across the entire overdoped region down to p 0.19 
where the pseudogap opens and the ratio then collapses 
rapidly, thus exposing an abrupt crossover to "weak" su- 
perconductivity as the Fermi surface reconstructs. 
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